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A COMPARISON BETWEEN OBSERVED WINDS AND CLOUD MOTIONS 
DERIVED FROM NIMBUS II HRIR MEASUREMENTS 

William E. Shenk 
and 

Earl R. Kreins, Major, USAF 
ABSTRACT 

Cloud motions over oceanic areas were obtained from HRIR measurements 
contained within the overlapped region of two adjacent Nimbus II orbits. These 
108- minute averaged motions were compared with instantaneous wind vectors 
interpolated from constant pressure charts. Corrections for uncertainties in 
the attitude of the satellite were made by superimposing coastlines that were 
within the fields of view for both orbits. Cloud mass centroids were determined 
to estimate the position for each cloud. Cloud top altitudes were obtained from 

t 

the Nimbus II HRIR data and the isotherms on the constant pressure charts. 
Generally, the direction and speed of the cloud motions were within 20° and 15 
knots of the observed wind direction and speed. 
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A COMPARISON BETWEEN OBSERVED WINDS AND CLOUD MOTIONS 
DERIVED FROM NIMBUS H HRIR MEASUREMENTS 
William E. Shenk 
and 

Earl R. Kreins, Major, USAF* 


BACKGROUND 

One of the primary objectives of the Global Atmospheric Research Program 
(GARP) is the determination of winds at several levels in the atmosphere and 
with 400 km horizontal spacing. Present and planned efforts to accomplish this 
task are by tracking constant pressure balloons and by measuring the speed and 
direction of cloud movements as detected by visible and infrared sensing devices 
on board satellites placed in geosynchronous orbit. To date, only visible sensors 
have been placed in synchronous orbit on board the Applications Technology 
Satellites (ATS), and the reported research results indicate that judiciously 
selected small clouds can serve as atmospheric tracers. The most widely em- 
ployed technique to obtain cloud motions from the ATS data is based on the ability 
to identify a cloud element on a series of successive photographs taken at inter- 
vals of approximately 20 minutes and then measure the cloud displacement from 
the first photograph to the Jast. Since the visible data provide no clear indication 
of cloud height, the clouds have been assigned to statistically determined preferred 

*Air Weather Service member temporarily attached to the Goddard Space Flight Center. Present 
affiliation: Joint Meteorological Satellite Program Office, Headquarters USAF, Washington, D.C. 
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pressure levels (one each for low and high clouds) or to the pressure level where 
the conventionally observed wind most closely approximates the cloud motion. 
Generally, the results from these techniques have been the determination of 
wind direction to within 20° and speeds to within 10 knots. Infrared sensors 
placed in geosynchronous orbit should provide more definitive information on 
cloud top heights. 

Most of these comparisons between cloud motions and wind observations 
have been performed within 50° of great circle arc of the subsatellite point. 
Figure 1 illustrates the geometric relations along a great circle passing through 
the subsatellite point between local zenith angle 8 , nadir angle 4> and great circle 
arc d. As the distance from the subsatellite point increases, a cloud is viewed 
more obliquely. As a result, greater difficulty is encountered in accurately 
identifying and locating a cloud element. Thus, it would be expected that cloud 
motions determined where 5 is large would be less accurate and compared less 
favorably with the observed winds than where 8 is small. 

Figures 2 and 3 show the portions of the Northern Hemisphere that would be 
covered by three and four equally spaced equatorial geosynchronous satellites, 
respectively, within three different limits of 8. Coverage over North America 
was purposely given low priority since this is one region where conventional 
meteorological data are adequate. With the system of three geosynchronous 
satellites, complete coverage at the equator is accomplished with 8 = 70°. If a 
four satellite system is employed, then complete coverage is assured to 30 8 N 
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within 8 = 60°. Let it be assumed that at 8 = 60° the capability to determine 
winds from the geosynchronous satellite has dropped to the point when alternate 
methods would receive strong consideration. Even with the four geosynchronous 
satellite coverage, cloud motion information would not provide wind estimates 
over large areas of the North Atlantic and North Pacific oceans. In the Southern 
Hemisphere, with its larger percentage of ocean surface to land surface, the 
oceanic data gaps would be more extensive. 

The concept of wind measurement from satellites was mentioned long before 
the launch of TIROS I (Widger, 1957). Greaves (1965) tried to implement the 
concept by recognizing that there is considerable overlap in the data from ad- 
jacent orbits at high latitudes for satellites with highly inclined orbits. Using 
TIROS IX video data, the cloud displacements from one orbit to the next were 
computed after corrections were applied for parallax and for the attitude differ - 
ences between orbits by matching recognizable landmarks in each orbit. These 
computed motion vectors were then compared with 500 mb geostrophic winds. 

The computed directions were generally within 30° of the observed direction ex- 
cept for three vectors where the average error was 140°. These three vectors 
were estimated from what were judged to be cirrus clouds well above the 500 mb 
surface. Generally, the computed speeds were substantially less than the 500 mb 
geostrophic winds with which they were compared. This effort suffered from the 
inability of the video data to locate the cloud at a specific level. 

The above discussion illustrates several key points. First, the visible data 
from the equatorial geosynchronous satellite may not be able to provide accurate 
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wind measurements at large local zenith angles and will provide only daytime 
information. Second, when the preferred levels technique is used, wind informa- 
tion is assigned to only two levels. Also, although the technique of obtaining 
winds from the video data in the overlapped region of the two adjacent orbits from 
a near polar orbiting satellite showed some promise, there was still an obvious 
need for better cloud top height estimates. Therefore, the ideal sensor to pro- 
vide winds at high latitudes would be a polar orbiting satellite with sufficient 
resolution to see some cloud detail and give information on cloud top height. The 
High Resolution Infrared Radiometer (HRIR) on board the Nimbus I, Nimbus II 
and Nimbus III satellites comes closest to this ideal. For this study, the Nim- 
bus II HRIR was chosen. Besides the search for a useful method for obtaining 
winds at high latitudes, this study can also examine the general question of the 
feasibility of obtaining winds from radiometric measurements. 

NIMBUS H SATE LUTE 

The Nimbus II satellite was launched into a near-circular retrograde sun- 
synchronous polar orbit with an inclination of 100.3°, apogee and perigee heights 
of 1179 km and 1095 km, respectively, and a nodal period of 108.17 minutes. A 
complete description of the spacecraft system is given in the Nimbus II Users' 
Guide (1966). 

HIGH RESOLUTION INFRARED RADIOMETER DATA 

The Nimbus II High Resolution Infrared Radiometer (HRIR) sensed upvelling 
radiation in the 3.5-4. 1 atmospheric "window" with a resolution w-hich was 8.6 by 


9.7 kilometers at the subsatellitc point from an orbital attitude of 1100 km. In 
this wavelength region, thermal emission and reflected solar radiation contribute 
about equally to the observed radiance during the daytime. Therefore, only at 
night is the detection of pure thermal radiation possible. Moonlight reflection 
from the earth is negligible. Thus, at night, reasonable meteorological interpre- 
tations of terrestrial radiation patterns can be made. The detected radiances 
are converted to equivalent blackbody temperatures. These derived tempera- 
tures are used to determine surface temperatures, under clear sky conditions, 
and cloud top heights where the clouds are assumed to be opaque and completely 
fill the view of the radiometer. 

The radiometer scanned perpendicular to the orbital track of the satellite. 
Space radiance levels are viewed along part of the scan while the radiometer is 
not pointing towards the earth. In the course of data processing and after con- 
verting the radiances to effective blackbody temperatures, the computer averages 
the effective blackbody temperatures over the space portion of the scan. All ef- 
fective blackbody temperatures that fall below this averaged value during the 
earth scan are given a value of 190 K. Since the average effective blackbody 
temperature associated with the space scans during the lifetime was 225 ± 5 K 
except for a short period foil' wing launch, any averaging of these temperatures 
with other effective blackbody temperatures above the space threshold might re- 
sult in fictiously low temperatures. Omitting the 190's would give a bias toward 
effective blackbody temperatures that are too warm. Thus, temperatures that 
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are near the threshold and where data from more than one scan spot were 
averaged must be cautiously interpreted. 

The infrared data currently are available in two forms. One is a photo- 
graphic image produced from the analog records. The other, a grid print map, 
is a computer display of the data rectified for various scales and map projections. 
More complete descriptions regarding these two types of data display can be 
found in the Nimbus II Users' Guide (1966) and in other publications such as 
Warnecke, et al. , (1968). 

The measurements from the IIRIR contain both periodic and random noise 
components. Spectral analysis revealed a strong 200 Hz noise component caused 
by interference from the spacecraft clock. This can be largely removed by 
numerical filtering, and a computer subroutine written by McMillin (1969) is 
patt of the standard mapping program. 

Williamson (1969) has investigated the random noise of the HRIR measure- 
ments. He calculated the noise equivalent temperature difference, NEAT, for the 
IIRIR at several temperatures. For a radiometric measurement, NEAT is de- 
fined as the change in target temperature which would produce a change in 
measured response equal to the r.m.s. noise at the output of the instrument. 
Figure 4 is reproduced from Williamson's paper and shows NEAT for a single 
digitized sample as a function of target temperature. It is apparent that for 
temperatures colder than 240 K the NEAT become large, even for data treated 
by McMillin's filter. In order to improve the temperature resolution at a given 
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grid point on a grid print map, the number of digitized samples, N, must be in- 
creased since the data error at any grid point is the error of any given sample 
divided by /N. Some data resolution must be sacrificed if improved temperature 
resolution is desired. 

The sun shade on the HRIR was not sufficiently large to prevent some addi- 
tional solar energy from affecting the measurements. These effects, as shown 
by the streaks in the pictorial presentation (Figure 5) between approximately 
55 N and 70 N, and in an analyzed grid print map (Figure G) for a small portion 
of Figure 5, are most noticeable when the sun was near the satellite horizon. No 
cloud motion calculations were performed within these regions. 

PROCEDURE 

The first step was to select two adjacent orbits with HRIR coverage at high 
latitudes so that there was considerable overlap in the two data fields. A Mer- 
cator map projection was used with a scale of 1:2 million. Thus, at the equator, 
this scale gives a distance of 15 nautical miles between grid points, while at a 
latitude of 60° the distance is approximately 8 nautical miles due to the geometry 
of the projection. Except at nadir angles very close to 50° at high latitudes, there 
was at least one measurement present at each grid point. Close to the subsatel- 
lite point, the reported effective blackbody temperature at each grid point was 
derived from the averaged data from several scan spots. 

From orbit to orbit there is a random error in locating the radiation data 
on the grid print maps. This error results from uncertainties in the attitude of 
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the satellite. These location errors can be as much as 30 nautical miles. This 
attitude error alone could result in sizeable differences between apparent cloud 
motions and the corresponding observed winds. A method tc reduce this error 
is to demand that a distinquishable land feature be detectable in both orbits. 

Then the coastlines, lakes, etc. can be matched and allowances made for the 
differences in attitude between orbits. 

Another location error is due to inaccurate determination of the relative 
location of successive scans (Sabatini and Sissala, 1969). This is caused by the 
errors in locating both horizons during the course of a scan. At the subsatellite 
point, this error is about 3 nautical miles and is somewhat greater nearer the 
horizons. Since this error is much smaller than the attitude error and since 
there is no simple technique of correcting for it over large areas, no correction 
was applied. 

Following the determination of the distance that one map had to be shifted 
to correct for attitude errors, individual cloud elements were selected that could 
be identified on both maps. Clouds that were in o:: very near surface frontal 
zones were avoided because of the expected sizeable influence of strong vertical 
motions near the front. Frontal zone clouds have the tendency to drift with the 
front and be poorly related to the winds moving through the front. Also, no 
cloud motion calculations were performed where the cloud top temperature was 
less than 240 K due to the nearness to the space threshold and the increase in 
random noise at cold temperatures. The centroid of each cloud mass was 
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computed and the movement of the centroid from the earlier to later map was 
called the computed cloud motion. 

The effects of parallax must be considered when an object that is closer to 
the viewing platform than the background is located on its background and is 
viewed from two different points in space. Figure 7 presents an exaggerated 
view. Since the clouds were geographically placed by referencing the earth's 
surface and the clouds are above the surface, the apparent zonal component of 
motion is more than the real zonal motion for the retrograde Nimbus II orbit. 

As the satellite moves through a full orbit, the earth rotates beneath it at a rate 
of 15° of longitude per hour. Therefore, at the equator, a point has moved about 
3010 km eastward during one Nimbus II orbital period. Thus, for a stationary 
cloud top at 6 km above the earth's surface with the satellite at an altitude of 
1100 km, there is an apparent 19 km easterly motion of the cloud. Thus, con- 
sidering the orbital period of the satellite, the apparent rate of movement would 
be 6 knots. At 60° N, this apparent rate of movement diminishes to 3 knots due 
to the smaller easterly shift of a point during an orbital period. Since most of 
the cloud motions computed were in the latitude range from 40° N - 60° N, no 
parallax correct n was applied unless a cloud top was placed above 500 mb. 
Then, a 5 knot westward component was added. 

There are several assumptions that are necessary in order to place a cloud 
top at a pressure level and to geographically place the centroid. It must be as- 
sumed that in the region where there is an apparent cloud that the cloud top 
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surface completely fills the field of view of the radiometer a"d it must be further 
assumed that the cloud is opaque. The failure of either assumption means that 
the apparent cloud top level is lower than the true level. To minimize the effect 
of these assumptions, the centroid was computed using the coldest region of the 
apparent cloud. A third assumption is that the cloud moves with the wind at the 
cloud top; in other words, the cloud is wafer thin. Serebreny et al. , (1969), com- 
paring cloud motions computed from ATS pictures and observed winds has indi- 
cated that this assumption is reasonable. 

Another consideration in the vertical placement of the cloud tops is the 
temperature corrections that could be applied to the HEIR data due to the inter- 
ference of atmospheric constituents. The principal interfering atmospheric 
constituents are carbon dioxide and water vapor. At the earth's surface over 
the range of nadir angles from 0° to 50°, corrections for these two constituents 
average about 3 K (Warnecke, McMillin and Allison, 1969). Naturally, above the 
surface where the masses of carbon dioxide and water vapor are less, the cor- 
rections would be smaller. Since 3 K corresponds to an average altitude error 
of about 0.3 km and most corrections would be smaller than 3 K, no corrections 
we” applied. 

Instead of assigning the coldest individual value at a grid point within the 
area used to compute the centroid, a spatial average of several points was taken 
as representative of the effective cloud top temperature. This further reduces 
the effect of random noise and gives ample weight to a sizable amount of the 
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cloud top surface. The pressure level of the cloud top was determined by match- 
ing the effective cloud top temperature to the same temperature determined by 
interpolating that temperature from constant pressure charts at the standard 
pressure levels. Then, the observed wind at the determined pressure level was 
interpolated from these conventional charts. 

Figure 8 shows the two data fields used to locate and identify a cloud ele- 
ment as well as the comparison between the computed cloud motion and the ob- 
served wind. 

DATA SELECTION 

In order to minimize the random errors in data location introduced by satel- 
lite attitude variations, a noticeable thermal difference between water and land 
was necessary because of the coastline matching that is essential to the pro- 
cedure. This nighttime thermal contrast between the land and the warmer water 
is greatest in the fall season during the six month period of available Nimbus II 
data. Therefore, all the cases selected were during the late summer and the 
fall seasons. Table 1 lists the orbital pairs, the midpoint time between the orbits, 
and the geographical features used to adjust for relative attitude errors. 

Synoptic situations were selected where most of the computed cloud motions 
were over water. Since the orbital period is 108 minutes, it becomes essential 
to minimize possible distortions in the cloud over this interval. Rough terrain 
influences changes in cloud shape, cloud height, and apparent cloud motions due 
to preferred regions of cloud development and dissipation. Not only do these 
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influences lead to erroneous cloud motions but they also increase the likelihood 
that a cloud will not be properly identified in the second orbit as the same cloud 
that was viewed in the first orbit. 

Over and east of Asia, the Daily Weather Map series supplied by the Japan 
Meteorological Agency (1966) was the conventional meteorological data source 
used to extract the wind measurements and the temperature information neces- 
sary to fix the pressure level of the cloud top. This series contained charts at 
the mandatory levels from surface to 500 mb. Above 500 mb, National Meteorologi- 
cal Center (NMC) charts supplied the same information. Over North America, all 
interpolations were made from NMC mandatory level charts. 

For the Asian synoptic cases, no observed wind interpolations were per- 
formed between the 1200 GMT constant pressure charts (closest time to the 
satellite passes) and the 0000 GMT constant pressure charts for the next day, 
except for one case. In this case, a rapidly moving short wave trough was prop- 
agating through the area causing large wind Auctions in a short time. There- 
fore, some adjustment was made to the 1200 GMT observed winds using the ob- 
served wind field displayed on the following set of 0000 GMT charts. The midpoint 
times between each orbital pair were generally around 1500 GMT (See Table 1). 

Space and time interpolation was performed for the North American cases 
since the midpoint time for the adjacent orbits was about 0600 GMT. 
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SAMPLE CASE 


A complete discussion of one of the seven cases will be presented followed 
by tabulation of pertinent data for the other six. The sample case to be discussed 
occurred November 11, 1966, with the midpoint time of 1530 GMT between the 
two adjacent orbits (data orbits 2400 and 2401). Figure 9 shows the amount of 
overlap in the HRIR data for the two orbits. Extensive portions of the China 
coastline and the Korean peninsula were clearly discernible in both orbits. A 
sharp 10 K temperature difference between land and water marked the position 
of the China coast and a lesser difference was observed to indicate the Korean 
coastline. These temperature gradients were then superimposed and the amount 
of attitude error between orbits determined. 

Table 2 lists the appropriate information of each comparison made between 
the seven computed cloud motions and observed winds. Figure 10 shows a 
schematic comparison of the computed vectors and the observed winds plus the 
determined cloud top levels. No adjustment in the observed winds was made be- 
tween the conventional map time (1200 GMT) and the midpoint time for the two 
orbits. Five of the seven computed cloud motions were from the simple trans- 
lation of a cloud element but the other two were attempts to use other techniques. 

The first of these was a very prominent cloud edge. Figure 11 shows two 
views of this cloud edge for the two orbits. Instead of a centroid computation, 
the position of the edge was indicated at the midpoint of the edge and in the 
strongest thermal gradient. Selection of the appropriate temperature to vertically 
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locate the cloud top was more subjective. This cloud edge appears to be the 
leading edge of a middle and high cloud shield moving off the Asian coast. Since 
cirriform clouds form at least part of the shield, it seems plausible that the 
radiometer is viewing a cloud that is not completely opaque. Therefore, the ob- 
served HPIR gradient from east to west is probably due to an emissivity gradient 
rather than a substantial change in the height of the cloud or the amount of cloud 
filling the field of view of the radiometer. Thus, more weight was given in the 
cloud height determination to the effective temperatures well behind the cloud 
edge. Table 2 indicates that 245 K was the average effective temperature selected. 
This estimate could easily be off by ±5 K. The interpolated cloud top level was 
450 mb and the computed cloud motion was 275 degrees at 65 knots with a 5 knot 
westerly component being subtracted due to parallax. The interpolated observed 
wind at 450 mb was 280 degrees at 70 knots. 

Another technique was the use of a hole in the clouds. This hole, shown in 
Figure 12 for the two orbits, was in the same general cloud mass that produced 
the prominent cloud edge to the northeast. The movement of the hole is con- 
trolled by the movement of the surrounding clouds. Therefore, the pressure 
"level” of the hole should be estimated from the surrounding cloud tops. A pres- 
sure level of 425 mb was assigned to the average effective temperature of 245 K 
and the cloud direction and interpolated wind direction agreed to within 15 de- 
grees while the cloud and observed wind speeds were the same (60 knots). 
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HIGHLIGHTS OF THE OTHER CASES 

Of the remaining six cases, four were near the Asian Coast and the other 
two over North America. All of the computed cloud motions in the Asian cases 
were over water and all of the cloud vectors computed over North America were 
over land. Table 3 give6 a complete list of the computed cloud motions for the 
six cases. In addition to the simple translation of cloud elements the centroid 
of an ensemble of small clouds was computed and the displacement of the en- 
semble measured. The comparison between the cloud motion and the concurrent 
wind was within 10 degrees and 10 knots. 

The greatest disparities between computed cloud motion direction and ob- 
served wind direction were noted in five computations which occurred over land 
between 850 mb and the surface. Four of the five cloud motions were computed 
for clouds over the Laurentian shield of eastern Canada and the fifth over the 
hilly terrain of northeastern Wisconsin. It is possible that there could be enough 
land influence on the clouds to cause the observed directional differences of 
30 - 40 degrees. These five vectors are marked with an asterisk in Table 3 and 
will be specially considered in the statistical analysis in the next section. 

STATISTICAL ANALYSIS 

Figure 13 shows the frequency distribution of the cloud top heights. The 
distribution appears to be trimodal, reflecting some tendency for clouds to be 
either low, middle, or high. 
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Figure 14 is a cumulative frequency diagram for the entire data sample 
which shows the numerical change in the effective blackbody temperature from 
the first to the second view of each cloud. There was only a 3 K change in the 
cloud top temperature in l r 8 minutes for over 60% of the clouds used as tracers. 
The majority of this amount of temperature fluctuation could be caused simply by 
viewing the cloud at different nadir angles from one orbit to the next. Over 85% 
of the clouds had an effective blackbody temperature change of <6 K. Of the five 
clouds with orbit to orbit effective blackbody temperature changes of >7 K, two 
of these were a cloud edge pair and the hole in the clouds. These results sug- 
gest that little convective activity was connected with the tracer clouds. This is 
not surprising since the clouds were sensed at high latitudes where minimal 
convective activity would be expected. The small temperature changes also 
indicate that non-con vective cloud top temperatures are conservative with time. 

Figures 15 and 16 respectively show plots of computed cloud direction 
versus observed wind direction and computed cloud speed versus observed wind 
speed. The 45° angle line in both figures represents a perfect correlation. The 
cloud motions where substantial topographical influence is suspected are labelled 
with an X. Since a high percentage of winds have a westerly component, there is 
a natural bunching of data points between 220° and 320° in the direction scatter 
diagram. The scatter appears to be about equally divided on either side of the 
perfect correlation line. Thus, no consistent directional bias in the computed 
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cloud motions is present. In the speed scatter diagram, there appears to be a 
slight tendency for the computed speeds to be higher than the observed speeds. 
This effect may be due to cloud propagation. More likely, it reflects the in- 
ability of the radiometer to verify that the field of view is completely cloud filled 
and that the cloud top is not completely opaque. These latter effects mean that 
some of the contributed radiance to the total sensed radiance is emitted from 
within or below the cloud. The results of these effects are that the true cloud 
tops are probably higher than those indicated by the radiometer. The opacity 
effect is probably greatest where the clouds are close to the cirrus level. Kulm 
and Weickmann (19G9) have reported that emissivities associated with cirrus 
are substantially less than 1.0. Thus, since the wind speeds normally increase 
with height below the tropopause, the observed wind speeds would have more 
closely approximated the cloud speeds at some higher level. 

Court (1958) and Durst (1954) have developed statistical relations to compute 
vector correlation coefficients. Court's total vector correlation coefficient, R wz 
is more complete. Let 

n 

Z> 

k* 1 

and 

n 

k 3 I 
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be two samples of winds where W = iU + jV and Z = iX + jY are wind vectors 
separated in time, in space, or in both time and space. The simple linear cor- 
relation coefficients, r ux , r yx , r uy , r yy and r are contained in the expres- 
sion for R WI where u, v, x and y arc deviations from the means U, V, X and 
Y. The sample variances of u and v are written as S u 2 and S v 2 , respectively. 
The expression for R wI 2 is given as 


S J « 'J, - 2r „, T „ r „ ) ^ S . J ( r , a , 

(S„ ,+ S,’)(«-r') 


The Durst coefficient, r wi 2 , is defined by 


r 

W I 



and the notation is the same that is used in the preceding paragraph. Court's 
relation involves multiple correlation of these variables while Durst's is like 
the simple correlation of two variables. Lamberth (1966) compared results of 
the two coefficients on the same data sample and found that the values of R wx 
are r x2 were similar except when R wi < 0.30. However, Lenhard (1967) dis- 
cusses some earlier work by Lenhard et al. , (1963) where R wi = 0.70 and 
r „ = 0.40. He cautioned against computing only r wi . 

Two computations of R^ 2 and r Wf were performed. One set of computations 
was with the full data sample (43 pairs) and the other excluded the five pairs 
where topography possibly influenced the clouds. The results for each sample 
were identical, R WI = 0.85 and r w2 = 0.84. The fact that R wx and r wi were the 
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same for the two samples shows that even though there was a relatively large 
error in direction between the low-level cloud motions and observed winds over 
hilly terrain, the small speed errors compensated. 

There is no rigorous way to estimate the confidence iimits for Court's R w . 
However, Court* suggested that a Fisher z transformation should be useful in 
testing the significance of linear correlation coefficients. In application of this 
test, comparison is then made with the normalized standard deviation, a , equal 
to 1/V2N-10, where N represents the number of independent data pairs. 

In the use of the z transformation, the null hypothesis is that the correlation 
between the two populations is zero. The probability that the observed correla- 
tion between the samples could be due to accidental sampling is determined. In 
this study, it is assumed that the distribution of sample correlation coefficients 
that could be computed by further studies is nol normal. This is a conservative 
view. The transformation, z, is given by z =1/2 [-i nil + r) - - r)] where, 

in this case s R W2 . As often is the case in meteorology, estimating the number 
of degrees of freedom is not straightforward. One sample contained 43 pairs 
and the other 38. However, the number of degrees of freedom is undoubtedly 
smaller due to spatial dependence between sample pairs within a given synoptic 
situation and some intercon elation between synoptic situations. The number of 
degrees of freedom as suggested by Court is 2N-10. A conservative hypothesis 
would be to assume that the seven synoptic situations are independent but that 

•personal communication, 1969 
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there is perfect correlation between pairs within each situation. A more reason- 
able degrees of freedom estimate probably lies in between this extreme and the 
extreme of complete pair independence. It was assumed that within any synoptic 
case, exactly two pairs were independent of one another (rather than the average 
number of six pairs). Thus, the number of the degrees of freedom, 2N-10, is 18. 
The student's t test relation is expressed as t = z/cr z which equals 5.48. For 
18 degrees of freedom this value of t indicates that the probability of the cor- 
relation coefficient originating from an uncorrelated population would be less 
than 0.1%. 

CONCLUSIONS 

The primary conclusion of this investigation is that by using Nimbus II HRIR 
data from two adjacent orbits, it is possible to obtain representative cloud mo- 
tions that correlate reasonably well with the observed winds. A substantial 
contribution to this positive result was undoubtedly the radiometer capability to 
more clearly determine cloud top height than is possible with visible data. This 
capability, at least in part, compensated for the poorer spatial and time resolu- 
tion than is currently available in computing cloud motions from the ATS satel- 
lite visible data. Since the overlap between adjacent orbits is most sig ificant 
at high latitudes, this technique might be most useful in the determination of 
winds over high latitude oceanic regions. It follows that it should be possible to 
derive wind data from other radiometric measurements where the resolution 
and noise characteristics are the same or better than the Nimbus II HRIR. This 
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is important since there are active programs to place radiometers in geosyn- 
chronous orbit with the wind measuring capability as the primary scientific ob- 
jective. Results as good or better than the above should be anticipated from 
geosynchronous altitude in view of the shorter interval between successive views 
of a cloud element. 
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Figure 2. Northern Hemispheric coverage with a system of three equatorial geosynchronous satel- 
lites. Curved lines are lines of equal 5 relative to each subsatellite point and the black dots in- 
dicate each subsateliite point. 
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Figure 3. Northern Hemispheric coverage with a system of four equatorial geosynchronous 
satellites. The nomenclature is the same as Figure 2. 
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Figure 4. Noise equivalent temperature difference (NEAT) for a single Nimbus II 

HRIR measurement. 
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Figure 7. Apparent cloud edge motion due to parallax for a satellite with a retrograde orbit. 
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Figure 10. Schematic presentation of the seven cloud motions and observed winds for the case 
of November 11, 1966. The cloud top heights are shown in mb. 








35 


Figure 12. Well defined hole in on extensive cloud sheet (Orbits 2400 and 2401). Isotherms 
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Figure 14. Cumulative frequency diagram of the orbit to orbit change in 
the effective blackbody temperatures of the cloud tops. 
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Figure 16. Scatter diagram of computed cloud motion speed versus observed wind speed. 
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TABLE 1 


List of orbit pairs, midpoint time between orbits, and geographic features 
used to correct for relative attitude errors 


Date 

Midpoint Time (GMT) 

Data Orbits 

Geographic Features 

September 15, 1966 

0530 

1636,1637 

Great Lakes 

October 15, 1966 

1425 

2040,2041 

Kamchatka Pen. 

October 28, 1966 

0535 

2208,2209 

Great Lakes 

November 1, 1966 

1545 

2267,2268 

China Coast 

November 8, 1966 

1330 

2360,2361 

China Coast 

November 10, 1966 

1410 

2386,2387 

Kamchatka Pen. 

November 11, 1966 

1530 

2400,2401 

Kamchatka Pen. 
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Pertinent data on cloud position 
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Pertinent data on cloud position, computed cloud top pressure level, 
computed and observed wind vectors, and cloud site and effective cloud 
top temperature on each arid print map for the remaining six cases 
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TABLE 3 - (CONTINUED) 
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